THE CARTAN GEOMETRY OF THE ROTATING KEPLER PROBLEM 
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Abstract. We investigate the Cartan and Finsler geometry of the rotating Kepler problem, a 
limit case of the restricted three body problem that arises if the mass of the one of the primaries 
goes to zero. We show that the Hamiltonian for the rotating Kepler problem can be regarded 
as the Legendre transform of a certain family of Finsler metrics on the two-sphere. For very 
negative energy levels, these Finsler metrics are close to the round metric, and the associated 
flag curvature is hence positive. On the other hand, we show that the flag curvature can become 
negative once the energy level becomes sufficiently high. 



1. Introduction 

The Kepler problem in a rotating frame arises as the limit of the restricted three body problem 
when the mass of one of the primaries goes to zero. Energy hypersurfaces of the Kepler problem 
as well as the restricted three body problem are noncompact. If one focuses on the bounded 
components of the energy hypersurfaces, i.e. one excludes the possibility of escape to infinity, then 
the only obstruction to compactness comes from collisions. However, it is well known in celestial 
mechanics that two body collisions can always be regularized. A particularly nice way to regularize 
the Kepler problem in an inertial frame was discovered by Moser in [lOj . By interchanging the roles 
of position and momenta Moser shows that the regularized planar Kepler problem is equivalent to 
the geodesic flow on the round two-sphere. Indeed, the point at infinity of the two-sphere corre- 
sponds to collision where the momenta explode. In pQ it was shown that for energies not too large 
the compact components of the regularized energy hypersurfaces for the planar restricted 3-body 
problem are starshaped. Here starshaped refers to the new cotangent bundle structure obtained 
after interchanging the roles of position and momenta; with respect to the old (g,p)-coordinates 
the Liouville vector field is given by q-^ and not by pj^. The challenging new question is now 
if starshaped can be replaced by the stronger condition of fiberwise convexity. If the energy hy- 
persurfaces are fiberwise convex we obtain a Finsler metric after Legendre transformation and the 
integral curves of the restricted 3-body problem can be interpreted as geodesies for the Finsler 
metric. The relevance of this problem comes from its relation with the project initiated by Helmut 
Hofer to construct finite energy foliations for the restricted 3-body problem. Indeed, if the energy 
hypersurfaces are fiberwise convex then the Conley-Zehnder indices of the closed characteristics 
correspond to the Morse indices of the corresponding geodesies and hence are nonnegative. Our 
first result is the following. 

Theorem A: The bounded components of the regularized planar Kepler problem in a rotating 
frame are fiberwise convex. 

In particular, the regularized Hamiltonian of the rotating Kepler problem is Legendre dual to 
a Finsler metric. The Legendre dual of a Finsler space was recognized by Miron as a Cartan 
space, see [11 J and the literature cited therein. Therefore, Theorem A provides the rotating Kepler 
problem with a Cartan structure. The sectional curvature of Riemannian geometry generalizes to 
the notion of flag curvature in Cartan, respectively Finsler geometry. In particular, for a Cartan 
surface the flag curvature is a function on the unit tangent bundle of the surface. This differs from 
the Riemannian case where the sectional curvature for a surface reduces to the Gaussian curvature 
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and hence is given by a function on the surface itself. Since in our case the underlying surface is a 
sphere, the flag curvature is a function on MP 3 regarded as the unit tangent bundle of the sphere. 
However, because the Kepler problem is invariant under rotations, the flag curvature is invariant 
under a circle action and quotienting out this circle action we again get a function on a sphere. 
However, be aware that this sphere does not coincide with the original one. The importance of the 
flag curvature comes from its relation with the Jacobi equation and hence with questions about 
indices, see [3]. Our next result is rather surprising since by Moser's result in an inertial frame, 
the regularized Kepler problem is equivalent to the geodesic flow of the round two sphere whose 
curvature is hence always positive. 

Theorem B: For the rotating Kepler problem the flag curvature can become negative. 

Let us explain our motivation for studying this problem. As mentioned above, the rotating Kepler 
problem is a limiting case of the restricted 3-body problem. Conley first proposed to use hyperbol- 
icity in the restricted 3-body problem in mission design, [HJ O |8] . Suppose we have the following 
problem. Given a satellite on an orbit around the earth, is it possible to move this satellite to 
an orbit around the moon using the gravitational forces of the earth, moon and maybe the sun 
without spending much fuel. It was Bclbruno who found the first realistic orbit using this idea, [J] . 
A variation of this orbit was succesfully put into practice to save the first Japanese lunar mission, 
0. 

Conley's idea takes advantage of a hyperbolic orbit around the first Lagrange point, also known 
as a Lyapunov orbit. These Lyapunov orbits only appear for sufficiently high energy. 

During a discussion with P. Albers, E. Belbruno, J. Fish and H. Hofer at the IAS in Princeton, the 
question arised whether (short) hyperbolic orbits can already arise below the first Lagrange point. 
For very small masses of the moon or for very negative energies it is known that this phenomenon 
never happens. Otherwise this question is completely open. 

A possibility to locate the position of hyperbolic orbits in phase space would be to find regions of 
negative curvature of the Finsler metric. Indeed, [S] contains a short computation which illustrates 
the relation between the flag curvature and the Conley-Zehnder index. In particular, one can 
clearly see that if the flag curvature is negative along an entire orbit, then this orbit is hyperbolic. 
We therefore hope that the methods studied in this paper have some applications for mission 
design. 

We conclude the introduction by pointing out that one could also try to understand the Conley- 
Zehnder indices by investigating whether there exists a convex embedding of a level set of the 
restricted 3-body problem into C 2 . The latter question has been answered in the affirmative in [5] 
provided the energy levels are sufficiently far below the first Lagrange value. 



2. A family of Minkowski metrics 
For p = (pi,jJ2) G M 2 and C > set p 1 - — (p 2 , — Pi) and consider the function H p : R 2 — — >• E, 

H p {q) :=(p ± ,q)-^+2C. 

Lemma 2.1. Suppose \p\ < C 2 . Then the curve H^ 1 ^) has two connected components: an 
unbounded one in the region {\q\ > 1/C}, and a bounded component T, p in the region {\q\ < 1/C} 
which bounds a strictly convex domain containing the origin. 

Proof. Given q ^ 0, let us determine all A € K — such that H p (\~ 1 q) = 0, i.e. 

\X\X-2C\q\X-\q\(p x ,q)=0. 
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Suppose without loss of generality that (p , q) > 0. Then the equation has the 3 solutions 



A» - C M + <g^j > o, A ± = -C|,| (i ± y'l - <^|j < 0. 

By hypothesis we have 

" |g|C 2 < ' 

and thus 1 - \A - x < 1 - (1 - x) < 1. We see that |A |, |A+| > C\q\ and |A_| < C|q|. Thus Aq x g 
and A^q belong to the bounded component contained in {\q\ < 1/C} and \Z 1 q to the unbounded 
component contained in {\q\ > 1/C}. 

For convexity we compute the gradient and Hessian 



Pi + |fpr y pw |g| 5 V -3qig 2 kh - 3q 



Vtfp(<7) = ^ J, lW/„(, y l = .-;-- 

A tangent vector to i?~ 1 (0) is given by 

v:=-(VH p (q)) ± = 
A short computation yields 

\q\ 5 (v,-RessH p (q)v) - (|q| 2 - 3q 2 ) (pi - - 6 qi q 2 (p 2 + T^s) (f i 

+ (M 2 -39 2 2 )(p 2 + ^) 2 

^H 2 M 4 + 2M(^,g) + l-3M 2 ( J5 ,g) 2 



To show positivity of the right hand side, we write (p^,q) = \p\\q\cost, (p,q) = \p\\q\ sin t. Then 
the right hand side equals \q\~ 2 f(t) with the periodic function 

f(t) := a 2 + 2acosi + 1 - 3a 2 sin 2 t, a := \p\\q\ 2 . 

The derivative f'(t) = — 2asint(l + 3acost) vanishes iff sint = or cost = — l/(3a). In the 
case sint = we have f(t) = (1 — a) 2 > 0, since < a < 1 due to the conditions \p\ < C 2 and 
\q\ < 1/C. In the case cost = — l/(3a) we have 

^ = ° ! -| + 1 - 3 ° 2 < 1 -(^> = 2 <5-° 2 >' 

The equation H p (q) = yields 

i - 2C = (p^,q) = \p\\q\ cost = - W = --L, 
M 3a 3\q\ 

thus |q| = 2/(3C). Using this and \p\ < C 2 , we estimate 

Hence a 2 < 1/3, which shows /(t) > and concludes the proof of the lemma. □ 

The convex curve E p endows T* M 2 with the structure of a Cartan space with fundamental function 
F* : R 2 ->• M> in the sense of [TT]: F*(q) is the unique A > such that \~ 1 q e S p . This A has 
been computed in the preceding proof to be 



F;(q) = C\q\ 1 + Wl 



kic 2 
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This Cartan structure is non-reversible (i.e. F*(—q) ^ F*(q)) for p ^ and does not appear to 
belong to one of the standard classes (Randers, Berwald, ...). 



3. Kepler problem in a rotating frame 

Consider the Kepler problem in a frame that is rotating around the origin with angular velocity 
a > 0. This system is again Hamiltonian with Hamiltonian 

1 \p\ 2 

H(p, q) = a{p 2 qi ~ Piq 2 ) - j^r + -y- 

= \{{Pi - Wif + (p 2 + a qi f) +U{q), 



with the potential 



The function U(r) = — - — ^-r 2 attains its maximum — |a 2 / 3 at the point r = a 2 / 3 . So the 
energy level 

{H(p,q) = -c}, c>^a 2 / 3 
has an unbounded component and a bounded component £ on which we have 

\q\ < a- 2 '\ 

The function 

H p (q)=H(p,q) + c 

is of the type considered above with p replaced by ap and 2C = |p| 2 /2 + c. Let us verify the 
hypothesis of Lemma |2.1[ 

Using c > |a 2 / 3 and setting x := \p\, it suffices to show 

Wax < (x 2 + 3a 2 / 3 ) 2 = x 4 + 6xV/ 3 + 9a 4 / 3 . 

Replacing x by xa 1 / 3 , this is equivalent to non- negativity of the function 

g(x) = x 4 + 6x 2 - 16x + 9. 

Since g attains its minimum at x = 1, the hypotheses of the lemma are satisfied. Hence the level 
set E is convex in q for each fixed p. Adding p = oo, E gives a fiberwise convex hypersurface in 
T*5 2 , where the coordinate on S 2 is p and the fiber coordinate is q. Let F* : T* S 2 — > K be the 
fundamental function of the corresponding Cartan space, which we computed above to be 



F*( P , q) = %f + 2c M 1 + Jl+ 16a{p± ^ 



4™ \ V \q\(\p\ 2 + 2c) 2 J ■ 

The fiberwise Legendre transform of F* (via ^F* 2 ) yields a Finsler metric F : TS 2 — > M whose 
geodesic flow equals the Hamiltonian flow of H on the energy level E. 

Note that H, and thus also F, is invariant under the rotation of S 2 around the vertical axis. 
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4. Polar coordinates 



To take advantage of the symmetry of F* it is useful to rewrite it in polar coordinates (a;, y) 
denned by 

Pi = x cos y, P2 = x sin y. 
The dual coordinates (r, t) on the cotangent bundle are determined by 

rdx + tdy — q\dp\ + q 2 dp2- 

In particular, 



/ x smy . cosy 

Ql = (cosy)r 1, q 2 = (smy)r H 1 

x x 



and hence 



1 



\q\ 2 =r 2 + ^t 2 , (p ± ,q) = -t 1 \p\ 



1 2 
= X . 



Hence in polar coordinates the fundamental function reads 
F*(x,y,r,t) = ]{x 2 + 2c)Jr 2 + -^t 2 (l + 1- 



IQat 



\ f 2 + ?' 2 (^ + 2c) 2 



Note that F* does not explicitly depend on y. 



4.1. A one-parameter family of Cartan metrics. The above formula for the Cartan metrics 
F* describes a two-parameter family with parameters c and a. However, each such Cartan metric 
is a constant multiple of one with a-parameter equal to 1. 



Proposition 4.1. F* a = a 1/3 F* a _ 2/3 r 



Proof. Use the coordinate transformation x — a 1 ' 3 ^', y — y' . On the cotangent fibers, this 
coordinate transformation induces r = a -1 / 3 / and t = t' . We get 



-< l 1 ^(x /2 + 2ca- 2 / 3 ) x /r> 2 



1 

z'V/ 3 



t ,2 \i+ 1 - 



IQat' 



\ y / r' 2 + ^i' 2 (a;' 2 + 2ca- 2 /3)2 / 



-« 1/3 C-/3, 1 (x',y',r',0- 
We shall henceforth put a = 1. 



□ 



5. Spray coefficients 

An elegant way of computing the flag curvature of a Finsler surface is via spray coefficients, see [31 
Section 12.5]. Since we are on the Cartan side we first have to Legendre transform the fundamental 
function F* , in order to apply the formulas in [3] directly. However, we do not know an explicit 
description of the Legendre transform of F* . We therefore derive in this section a formula to 
obtain the spray coefficients of a Cartan surface directly from the fundamental function. This 
allows us in turn to implement a MAPLE algorithm for computing the flag curvature of a Cartan 
surface. We abbreviate 

C* 



and denote by C the Legendre dual of C* . The coordinates on the tangent bundle are given by 

u(x,y,r,t) = C*(x,y,r,t) 
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v(x,y,r,t) = C* t {x,y,r,t). 

Here a subscript denotes differentiation with respect to the corresponding variable. We first recall 
some useful facts of the Legendre transformation for which we refer to pH Section 14.8] or [TT] . 
Since C* and C are 2-homogeneous they numerically coincide 

C*(x,y,r,t\ = c(x,y,u(x,y,r,t),v(x,y,r,t)y (1) 

The cometric coefficients are defined by 

g n =£* rr , g 12 =g 21 =Ct, g 22 = C* tt . (2) 

The metric coefficients inverse to the cometric coefficients can be obtained by 

.911 = £uui g\2 — 321 = £mij 322 = £vv 

Moreover, since the Legendre transform is involutive, one has the identities 

f £'Uj t by 

We further abbreviate the determinant of the cometric by 

g = g 11 g 22 -(g 12 ) 2 - 

Differentiating ([I]) we obtain 

dC* _ dC _ dC + dCdu dC dv 
dx dx dx du dx dv dx 



and therefore 

Similarly 

From the identity 

one obtains 



9 

In a similar vein, one derives the formulas 



c — C* — C C* — c c* 

'-x — >-x '-u'-rx '-v'-'tx- 
y ^u^ry ^"U^ty' 

Q _ <9r _ dCu dCudu dC u dv 
dx dx du dx dv dx 



-C C* — C C* 

-giiC* x - 3i2^L 

g C rx + g C tx 



— —C C* — C C* 



-(- 9 22 c; y + g 12 c* ty 



-C C* — c c 



tx 



r — —C r* — r r* 

= l -U 2 c: v -g^c ty 
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Plugging these identities into the formulas for the spray coefficients G, H (this H is not the 
Hamiltonian!) in [3, p 330] we get 

2G g £ vv £ x g/l vx /l v gH uv /l y -\~ gfL uy [L v 

~ 9 \J"x ~~ £u£>rx — f-'vC-tx^j — ^rx — .9 ^tx J 

+g 12 (c; - c u c; y - c v c* ty ) + ( - g 22 c* ry + g 12 c* ty )c v 

= (g^Cl+g 12 ^) - (g^C^+g 12 ^)^ 

, f 11 /»* _ 12 p* _ 11 /»* _ 12 22 /»* , 12 /»* \ /» 

^tx 9 ^rx 9 ^tx 9 t-'ty 9 L'ry .9 ^iy J 

= [ g n c x + 3 12 /:;) - (g ll c* TX + g 12 c* ry )c u - (g 12 c* rx + g 22 c; y )c v 



9 ^uu^y 9^uy£*u 9^vu^x 9^"UX^"U 

+9 f Ac — £"u,£ r x ~ £v£-tx^j + (.9 Ars — 5 ^tx^^u 

= (g 12 C* x + g 22 C* y ) - (g 12 C* tx + .g 22 /^) £„ 

i/' 22/.* _ 12 />* 22/.* _„12/.* , 12 f* _ 11 /.* \ /• 

T ly *-rj/ 5 1 *-tj/ y *-rs/ y L 'rx<9 J -rx 9 J "tx) L "u, 

= (g 12 C x + ,g 22 A*) - (g 12 C tx + g 22 C ty ) C v ~ (g 1 1 C* tx + g 12 C* ty ) £ 



5.f. A formula for the flag curvature following Bao-Chern-Shen. Wc shall use formula 
(f2.5.14) from [5] to compute the flag curvature. This formula is given by 

_ (G xv — G yu )v + 2GG UU + 2HG UV — G U G U — G V H U 

VCy 

All these terms can be rewritten entirely into terms that exist on the Cartan side. In particular, 
we shall use the above formulas for G and H in terms of g l i and C* and their derivatives. 

As before, we define the Hamiltonian or Legendre dual of C by 

1 9 

The cometric can be obtained via Equation As the Hamiltonian is given in terms of the base 
coordinates x, y and cotangent fiber coordinates r, t, we need to apply the chain rule to compute 
([3]), which is given in terms of x,y and the tangent fiber coordinates u, v. We need to use the 
following relations, 

d dr d dt d d d 

Similarly, 



du du dr du dt ^ U dr ® 21 dt 



d d d 

•5- = .912 -W- + .922 ^7 ■ 

dv dr dt 

Note that the metric can be expressed in terms of the cometric and its determinant. The depen- 
dence on the base coordinate is a little tricky, since the independent variables are now x,y,u,v, 
whereas the Hamiltonian is given in terms x, y, r, t. This means that the (co)-metric comes in. Let 
us write out the required expression. We shall use the notation 

dG 
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to indicate the derivative of G keeping y, u, v fixed. We obtain 

<'■'■)„ .,,- \ dx )y,r,t \ dr Jx,y,t\ dx )y,u,v V & ) X ,y,r\ dx 

dx )y,r,t \ dr Jx,y,t \\ dx Jy,u,v V dx ) y,u,v , 

9t )*,V,r (( dx \,u,v U+ ( dx \,u,v") ' 

Similarly, we have 

(K\ = (OG^ + (8G} ((Bn\ u+ fd 9l 



\ J x,u,v V J x.r,t V ) x,y,t V ^ } x,u,v \ / x,u,v 

8G\ ffdg 21 \ ^(dg 22 \ \ 

/ x,y,r \ \ & / x,u,v \ y / x, it, v 1 



In the above we have used 



5n u + QV2V an d t = g 2 \u + g 22 v. 



Using the cometric g tJ , we can conversely express u and v in terms of r and t. Hence we can 
translate all derivatives in the expression for the flag curvature in terms of quantities on the 
Cartan side. 



6. Computing the Flag curvature of the rotating Kepler problem 

We now come to the computer implementation of the curvature computation. The MAPLE- 
program is included in the appendix. Although the program itself is straightforward with the 
above formulas, there are some practical issues involved which we shall briefly discuss. 

Even though our computer setup was sufficient to compute the full expression for the flag curvature 
of the Cartan metric F*-^ as a function of x, r, t, the resulting expression was too unwieldy to 
simplify to reasonable proportions. However, it is possible to evaluate and simplify the curvature 
exactly in specific points, for instance if the r coordinate is 0. The simplified expression for the 
curvature along this ray is still very long, see Equation Q. 

In order to make these evaluations work, one only needs to realize that it is beneficial for a 
computer to quickly evaluate the results as soon as further differentiation is not needed. This is 
the reason for the many substitutions in the MAPLE program. 



6.1. Analytical and numerical results. Wc first fix the fiber coordinate (r,t) — (0,x) and 
compute the curvature as a function of x, which is the coordinate on the base sphere S* 2 . The 
MAPLE program yields the following expression for the flag curvature. Write 

a := V 'x 4 +4x 2 c + Ac 2 - 16 x 
to simplify the resulting expression. 

K(F* 1 )(x, 0, 0, x) = r 

(x 1 + 2 c + a) (x 2 a + 2 ca + x 4 + 4 x 2 c + 4 c 2 - 8 x) (x 4 + 4 x 2 c + 4 c? - 16 x) 

(4) 
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(5824 x 2 c 4 - 5888 x 3 c 5 - 3840 x 2 c - 2240 c 6 x - 6320 x 5 c 4 - 384 x 2 a + 1120 x 6 c 5 + 2 x u c 

+28 x 12 c 2 - 6528 x 5 c + 256 c 8 - 864 ac 5 x - 1872 acV + 896 x 2 c 7 - 1296 x 7 + 204 x 10 

-768 c 5 - 9 x 13 + 2096 x 8 c - 160 x 14 c - 1060 cV - 3520 cV + 11520 x 4 c 3 + 3840 c 3 a; 

+7584 x 6 c 2 - 5952 x 3 c 2 - 648 x 7 ac 2 - 126 x 9 ac - 1120 i 3 qc + 2448 x 4 ac 2 + 1152 ac 2 a; 

+1920 x A + 168 x 10 c 3 + 1344 x 4 c 6 + 560 x 8 c 4 + 1032 s 6 ac - 1584 x 5 ac 3 + 1632 ac 3 x 2 

+128 ac 7 + 384 ac 6 a; 2 - 9 x xl a + 2 x 12 ac + 132 x 8 o: + 320 x 6 ac 4 + 120 i 8 ac 3 

+24 x 10 ac 2 + 480 ac 5 ! 4 - 528 ax 5 - 384 ac 4 ) 

Since this expression is too complicated to deduce anything meaningful by hand, we have plotted 
some graphs in Figures [lj [2] and [3] We see that if c is small enough, then the flag curvature can 
become negative. Note that the graphs are also plotted for x < 0, where the sign of t changes 
since we insert the covector (r, t) — (0, x). 



in ; 

2U : 



-10 



0: 
-20 : 
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Figure 1. Flag curvature as a function of x for r = 0, t — ±1 and c = 1.51 




Figure 2. Flag curvature as a function of x for r = 0, t = ±1 and c = 1.65 




Figure 3. Flag curvature as a function of x for r = 0, t = ±1 and c = 2 

Let us now consider the flag curvature of K c as a function on the unit cotangent bundle. We shall 
use the coordinates x, cj>, where the cotangent fiber coordinates are given by 

(r, i) = (sine/), cos 4>). 

The full expression for the curvature for r ^= is considerably more complicated than Equation Q, 
so we shall only present numerical results for this general case. 

Using MAPLE we made a graph of the flag curvature by evaluating the curvature on a 256 x 256 
lattice, see Figure [4] On this lattice the minimal value of the flag curvature for c = 1.55 was 
—5.55, and the maximal value 15.21. 
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Figure 4. Flag curvature as a function of x and <fi for c = f .55 



7. Appendix 

Here is a general program to (numerically) compute the flag curvature of a Cartan metric. 

> K:=proc(F,x0,y0,r0,t0) 

> local L,Lx,Ly,Lr,Lt,Lxt,Lyr,Lrr,Lrt,Ltt,D, 

> G , H , Gr , Gt , Gar , Gat , Grr , Grt , Gxt , Gyr , Hr , K , 

> Garre , Garte , Grre , Grte , Gxte , Gyre , Hre , Gre , Gte , 

> Ht,Hrr,Hrt,Htt,detg,gi,g,gie,ge, 

> Lsx , Lsy , Lsqlx , Lsqly , Lsq2x , Lsq2y : 
> 

> ## Hamiltonian given 

> H:=F~2/2: 
> 

> Hr:=diff (H,r) : 

> Ht:=diff (H,t) : 

> Hrr:=diff (Hr,r) : 

> Hrt:=diff (Hr,t) : 

> Htt:=diff (Ht,t) : 
> 

> gi[l,l] :=Hrr: gi [1 , 2] : =Hrt : gi [2, 1] :=Hrt : gi [2 , 2] : =Htt : 

> detg : =gi [1 , 1] *gi [2 , 2] -gi [1 , 2] '2 : 
> 

> ##obtain Finsler 'metric' g_{ij}: note that the dependence on the fiber 

> ##coordinates is not correct, since we haven't performed a Legendre 

> ##transf ormation 

> g[l,l] :=l/detg*gi[2,2] : g[l , 2] : =-l/detg*gi [1 ,2] : 

> g[2,l] :=-l/detg*gi[2,l] : g[2,2] :=l/detg*gi [1, 1] : 
> 

> Lr:=r: 

> Lt:=t: 
> 

> Lsx:=diff (H,x) : Lsy :=diff (H,y) : 
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> Lsqlx:=diff (Hr,x) : Lsqly :=diff (Hr,y) : 

> Lsq2x : =dif f (Ht , x) : Lsq2y : =dif f (Ht , y ) : 
> 

> ##Formulas from Section 5 

> G:= 1/2* ( gi[l,l]*Lsx+gi[l,2]*Lsy-Lr*(gi[l,l]*Lsqlx+gi[l,2]*Lsqly) 

> -Lt*(gi[l,2]*Lsqlx+gi[2,2]*Lsqly) ): 

> H : =1/2* (gi [1 , 2] *Lsx+gi [2 , 2] *Lsy-Lt* (gi [1 , 2] *Lsq2x+gi [2 , 2] *Lsq2y) 

> -Lr*(gi[l,l]*Lsq2x+gi[l,2]*Lsq2y) ): 
> 

> ##avoid double computations 

> Gar:=diff (G,r) : Gat :=diff (G,t) : 
> 

> Gr : =Gar*g [1,1] +Gat*g [2,1]: 

> Gt : =Gar*g [1,2] +Gat*g [2,2]: 
> 

> ##evaluate quickly 
> 

> ge [1 , 1] : =subs (x=xO , y=yO , r=rO , t=tO , g [1 , 1] ) : 

> ge [1 , 2] : =subs (x=xO , y=yO , r=rO , t=tO , g [1 , 2] ) : 

> ge [2 , 1] : =subs (x=xO , y=yO , r=rO , t=tO , g [2 , 1] ) : 

> ge [2 , 2] : =subs (x=xO , y=yO , r=rO , t=tO , g [2 , 2] ) : 

> gie[l,l] :=subs(x=xO,y=yO,r=rO,t=tO,gi[l,l]) : 

> gie[l,2] :=subs(x=x0,y=y0,r=r0,t=t0,gi[l,2]) : 

> gie[2,l] :=subs(x=x0,y=y0,r=r0,t=t0,gi[2,l]) : 

> gie [2,2]: =subs (x=xO , y=yO , r=rO , t=tO , gi [2 , 2] ) : 
> 

> ##avoid double computations 

> Garre:=subs(x=xO,y=yO,r=rO,t=tO,diff (Gr,r) ): 

> Garte:=subs(x=xO,y=yO,r=rO,t=tO, diff(Gr,t) ): 
> 

> ##evaluate when no more dif f erentation is needed. 

> Grre:=Garre*ge[l,l]+Garte*ge[2,l] : 

> Grte:=Garre*ge[l,2]+Garte*ge[2,2] : 
> 

> 

> Gxte:=subs(x=xO,y=yO,r=rO,t=tO, diff(Gt.x) ) 



> +subs(x=xO,y=yO,r=rO,t=tO, diff(Gt,r) ) 

> *subs(x=xO,y=yO,r=rO,t=tO, ( dif f (g [1 , 1] ,x) 

> *(gie[l,l]*r+gie[l,2]*t)+diff (g[l,2] ,x)*(gie [2, 1] *r+gie [2,2] *t) ) ) 

> +subs(x=xO,y=yO,r=rO,t=tO, diff(Gt,t) ) 

> *subs(x=xO,y=yO,r=rO,t=tO, ( dif f (g [2 , 1] ,x) 

> *(gie[l,l]*r+gie[l,2]*t)+diff (g[2,2] ,x)*(gie [2, 1] *r+gie [2,2] *t) ) ): 
> 

> 

> Gyre :=subs(x=xO,y=yO,r=rO,t=tO, dif f(Gr,y) ) 

> +subs(x=xO,y=yO,r=rO,t=tO,diff (Gr,r) ) 

> *subs(x=xO,y=yO,r=rO,t=tO, ( dif f (g [1 , 1] ,y) 

> *(gie[l,l]*r+gie[l,2]*t)+diff (g[l,2] ,y)*(gie [2, 1] *r+gie [2,2] *t) ) ) 

> +subs(x=xO,y=yO,r=rO,t=tO,diff (Gr,t) ) 

> *subs(x=xO,y=yO,r=rO,t=tO, ( dif f (g [2 , 1] ,y) 

> *(gie[l,l]*r+gie[l,2]*t)+diff (g[2,2] ,y)*(gie [2, 1] *r+gie [2,2] *t) ) ): 
> 

> 

> Hre:=subs(x=xO,y=yO,r=rO,t=tO,diff (H,r) )*ge[l,l] 
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> +subs(x=xO,y=yO,r=rO,t=tO,diff (H,t) )*ge[2,l] : 
> 

> 

> Gre:=subs(x=xO,y=yO,r=rO,t=tO,Gr) : 

> Gte:=subs(x=xO,y=yO,r=rO,t=tO,Gt) : 
> 

> K:= subs(x=xO,y=yO,r=rO,t=tO, 

> ( (Gxte-Gyre)*(gie[2,l]*r+gie[2,2]*t)+2*G*Grre+2*H*Grte-Gre*Gre-Gte*Hre ) 

> /( (gie[2,l]*r+gie[2,2]*t)*Lt) ): 
> 

> 

> end proc; 

> a:=l:c:=1.55:x0:=0.75: y0:=0: r0:=0: tO:=l: 

> F:=l/4*(x~2+2*c)*sqrt(r~2+x~(-2)*t~2)* 

( l+l*sqrt(l-16*a*(t)/sqrt(r~2+x~(-2)*t~2)/(x~2+2*c)~2 ) ); 

> K(F,xO,yO,rO,tO) ; 

The following program is optimized for the rotationally symmetric Cartan metric associated with 
the rotating Kepler Hamiltonian. Here we use the same notation as in Section [4] 

> K:=proc(F) 

> local H,Hr,Ht,Hrr,Hrt,Htt,gi,detg,g, 



> Lr,Lt,Lsx,Lsqlx,Lsq2x, 

> G , Gr , Gt , Gar , Gat , Grre , Gtt , Gre , Gt e , Hre , xO , rO , tO , 

> Garre,Garte,Grte,Gxte, 

> gie,detge,ge,Lsqlxe,Lsq2xe: 
> 



> ## Hamiltonian given 

> H:=F~2/2: 
> 

> Hr:=diff (H,r) : 

> Ht:=diff (H,t) : 

> Hrr:=diff (Hr,r) : 

> Hrt:=diff (Hr,t) : 

> Htt:=diff (Ht,t) : 
> 

> gi[l,l] :=Hrr: gi [1 , 2] : =Hrt : gi [2, 1] :=Hrt : gi [2 , 2] : =Htt : 

> 

> detg : =gi [1 , 1] *gi [2 , 2] -gi [1 , 2] "2 : 

> 

> ##obtain Finsler 'metric' g_{ij}: note that the dependence on the fiber 

> ##coordinates is not correct, since we haven't performed a Legendre 

> ##transf ormation 

> g[l,l] :=l/detg*gi[2,2] : g[l , 2] : =-l/detg*gi [1 ,2] : 

> g[2,l] :=-l/detg*gi[2,l] : g [2, 2] : =l/detg*gi [1 , 1] : 
> 

> 

> Lr:=r: 

> Lt:=t: 
> 

> Lsx:=diff (H,x) : 

> Lsqlx:=diff (Hr,x) : 
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> Lsq2x:=diff (Ht,x) : 
> 

> ##Formulas from Section 5 

> G:= 1/2* ( gi[l,l]*Lsx-Lr*(gi[l,l]*Lsqlx)-Lt*(gi[l,2]*Lsqlx) ): 

> H:=l/2*(gi[l,2]*Lsx-Lt*(gi[l,2]*Lsq2x)-Lr*(gi[l,l]*Lsq2x) ): 
> 

> ##avoid double computations 

> Gar:=diff (G,r) : Gat :=diff (G,t) : 
> 

> Gr : =Gar*g [1,1] +Gat*g [2,1]: 

> Gt : =Gar*g [1,2] +Gat*g [2,2]: 
> 

> ##evaluate quickly 

> xO:=x: r0:=0: tO:=x: 

> ge [1 , 1] : =subs (x=xO , r=rO , t=tO , g [1 , 1] ) : ge [1 , 2] : =subs (x=xO , r=rO , t=tO , g [1 , 2] ) : 

> ge [2 , 1] : =subs (x=xO , r=rO , t=tO , g [2 , 1] ) : ge [2 , 2] : =subs (x=xO , r=rO , t=tO , g [2 , 2] ) : 

> gie[l,l] :=subs(x=xO,r=rO,t=tO,gi[l,l]) : gie[l,2] : =subs(x=xO ,r=rO ,t=tO ,gi [1 ,2] ) : 

> gie[2,l] :=subs(x=xO,r=rO,t=tO,gi [2, 1] ) : gie[2,2] :=subs(x=xO,r=rO,t=tO,gi [2,2] ) : 

> detge : =subs (x=xO , r=rO , t=tO , detg) : 

> Lsqlxe : =subs (x=xO , r=rO , t=tO , Lsqlx) : 

> Lsq2xe : =subs (x=xO , r=rO , t=tO , Lsq2x) : 
> 

> ##avoid double computations 

> Garre:=subs(x=xO,r=rO,t=tO,diff (Gr,r) ): 

> Garte:=subs(x=xO,r=rO,t=tO, diff(Gr,t) ): 
> 

> ##evaluate when no more dif f erentation is needed. 

> Grre:=Garre*ge[l,l]+Garte*ge[2,l] : 

> Grte:=Garre*ge[l,2]+Garte*ge[2,2] : 
> 

> Gxte:=subs(x=xO,r=rO,t=tO, diff(Gt.x) ) 

> +subs(x=xO,r=rO,t=tO, diff(Gt,r) ) * (-ge [1 , 1] *Lsqlxe-ge [1 ,2] *Lsq2xe) 

> +subs(x=xO,r=rO,t=tO, diff(Gt,t) ) * (-ge [1 ,2] *Lsqlxe-ge [2 ,2] *Lsq2xe) : 
> 

> Hre:=subs(x=xO,r=rO,t=tO,diff (H,r) )*ge[l,l] 

> +subs(x=xO,r=rO,t=tO,diff (H,t) )*ge[2,l] : 
> 

> Gre:=subs(x=xO,r=rO,t=tO,Gr) : 

> Gte:=subs(x=xO,r=rO,t=tO,Gt) : 
> 

> K:= subs(x=xO,r=rO,t=tO, 

> ( (Gxte)*(gie[2,l]*r+gie[2,2]*t)+2*G*Grre+2*H*Grte-Gre*Gre-Gte*Hre ) 

> /( (gie[2,l]*r+gie[2,2]*t)*Lt) ): 
> 

> end proc ; 

> F:=l/4*(x~2+2*c)*sqrt(r~2+x~(-2)*t~2)* 

( l+l*sqrt(l-16*a*(t)/sqrt(r~2+x~(-2)*t~2)/(x~2+2*c)~2 ) ); 

> GK:=K(F) ; 

> curv:=simplify(GK,assume=positive) ; 

> plot({subs(a=l,c=2,curv)},x=-10. . 10) ; 
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